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Abstract 

We present recent results for heavy flavor and jets produced in ^snn = 200 GeV p-p, d-Au, 
Cu-Cu and Au-Au collisions at at RHIC. We find J/i/) production is supjircsscd in Au-Au, but 
high energy J/-i/) are not suppressed in d-Au or Cu-Cu. Non-photonic electrons from D and B 
mesons are both suppressed in Au-Au. We study jets using two-particle 2D angular correlations 
as well as jet reconstruction algorithms. Wc find that jets show little if any suppression, but 
are highly modified in central Au-Au, becoming elongated in the rj direction and having fewer 
high-pt particles but many more low-pt particles for central Au-Au. 
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1. Introduction 

The RHIC is a versatile accelerator that has been used to study p-p, d-Au, Cu-Cu 
and Au-Au collisions. In the present paper we are interested in the system produced 
in y^NN = 200 GeV Au-Au collisions. Heavy flavor and jets are produced early in the 
collision and we measure how their properties are modified by the medium produced in 
Au-Au collisions, using p-p and d-Au collisions to determine the baseline expectation. 

The primary component of the STAR detector is a TPC measuring charged tracks 
with pt > 0.15 GeV/c, covering 27r in azimuth and — 1 < < 1 in pseudorapidity. 
At larger radius is a barrel electromagnetic calorimeter, containing a shower maximum 
detector, which is used to measure neutral energy and identify electrons. The Forward 
Meson Spectrometer (FMS) measures 7r°s for 3 < r?^o < 4. The data acquisition system 
has been significantly enhanced, allowing data rates of up to 1000 events per second. The 
Time of Flight (TOP) system installation was completed before the most recent run. 
The TOF covers the same angular range as the TPC and increases the identification of 
7r^ and K='= to about 1.6 GeV/c and protons to 3 GeV/c as well as identifying electrons 
where the dE/dx crosses the tt* and bands near 0.3 and 0.7 GeV/c. An example of 
an analysis requiring the TOF is shown in the left panel of Fig. 1 
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2. Heavy Flavor 

Charm and bottom quarks are produced early in the coUision. Most charm and bottom 
quarks emerge as D and B mesons, so-called open heavy flavor, which have cr values of 
hundreds of microns. Wc would like to determine the c and b quark transport properties 
in the medium by measuring their energy loss and flow. Some charm and bottom quarks 
are bound into J/^) and T and their excited states, so-called quarkonia. It has been 
predicted that quarkonia will "melt" in a quark gluon plasma [1], the order of melting 
depending on the binding energies and plasma temperature. 



Quarkonia states with the lowest binding energy, the ■i/'(2S) and T(3S), melt at a 
lower temperature than the ,J/V; and T(2S) which melt at a lower temperature than 
the T(1S) [2]. At RHIC the T(1S) should not melt while the J/V' and T(2S) should be 
partially melted and the T(3S) should be almost completely melted. This melting would 
result in fewer observed quarkonia than expected from extrapolation of p-p collisions. 

We observe quarkonia states through the e+e" decay channel, identifying high energy 
electrons by using p/E in the barrel electromagnetic calorimeter and shower max detector. 
We measure suppression in A- A collisions using Raa = crAA/NbinaryCpp, assuming the 
production mechanism depends on hard scattering processes. Normally this is measured 
as a function of pt (so a shift in the spectrum would be interpreted as a suppression). 
For d-Au we see no strong evidence for suppression of i/tp. For Cu-Cu STAR measures 
the Z/tp cross section above 5 GeV and sees no suppression [3]. PHENIX does see a 
suppression at lower energies [4]. For Au-Au we observe a suppression with Raa ~ 
0.5 for i/tj} momenta below 5 GeV/c. 

The T is harder to measure and in the e+e~ channel it is diflicult to distinguish the IS, 
2S and 3S states so wc; combine them in our measurement. We have significant signal for 
p-p, d-Au and Au-Au collisions at = 200 GeV, finding R^au = 0.78±0.28±0.2. For 

the most central 60% of Au-Au collisions we measure Raa = 0.78 ±0.32 ±0.22. The large 
dataset taken during the last year with less material inside the TPC and including the 
TOF will reduce the uncertainties on Raa and we may be able to tell if the suppression 
is consistent with the 2S and 3S melting while the IS stays intact. 

Raa may change for reasons other than quarkonia melting. It is possible recombination 
is important and/or the feed-down to J/ijj from B meson decay could change. For p-p 
collisions we measure this feed-down using angular correlations of hadrons with respect to 
the J ftp direction [3]. Using Pythia and identifying J/ips that came from B meson decay 
versus directly produced we get two histograms. For B meson decays there is a same-side 
peak which is absent for directly produced J/^ps. We adjust the relative weights of these 
histograms to match the data and determine that the feed-down fraction, shown in Fig. 1, 
varies from about 10% to over 20% depending on J/ ip momentum. This is consistent with 
much higher beam energies, suggesting that once a threshold is crossed for production of 
B mesons then the ratio of b to c production is constant. If this feed-down were modified 
in A- A collisions Raa of the J/ip would be affected. 



2.1. J/tp, T at 




= 200 GeV 
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Fig. 1. Left panel shows the di-electron invariant mass spectrum from 200 GeV p-p collisions. The peaks 
are the oj, <p and 3/4>- The solid line is a cocktail of known contributions to the Mj,+^- spectrum which 
describes p-p well. The middle panel shows the fraction of Jip that come from decay of B mesons, 
comparing p-p collisions at -ysjviv = 200 GeV to higher energy collisions [5] . The right panel shows the 
fraction of non-photonic electrons from B meson decay for p-p collisions at ^sjvjv = 200 GeV. 

2.2. Open Heavy Flavor 

The cleanest way to measure open-heavy-flavor production would be to reconstruct 
D and B mesons, although this is difficult due to small branching ratios of measurable 
channels. We use the semi-leptonic decays, measuring single electrons. Most high energy 
electrons observed in heavy ion collisions come from 7 conversion to e~''e~. We subtract 
those using the small opening angle of the 7 e+e~, {6 m^/E^). After subtraction 
we are left with the non-photonic electrons, a mixture of D and B meson decays. Using 
Pythia we calculate electron-hadron correlations for D and B decays (the heavier B 
has a flatter near side correlation) and adjust the relative amplitudes to measure rb = 
NeB/(NeD + Neg) as a function of electron momentum. Above pt = 5 GeV/c we find 
rb « 0.5 as shown in Fig. 1 [6]. It has been observed that Raa for non-photonic electrons 
in ^snn = 200 GeV Au-Au collisions is suppressed to about the same level as inclusive 
hadrons [7], suggesting that charm and bottom quarks lose energy. It had been believed 
that heavy quarks would have little radiative energy loss and would not be suppressed [8] . 
The suppression of non-photonic Raa hi Au-Au collisions combined with rb measured in 
p-p collisions means that Raa ^-^ ^'^'^ confidence level. 



2.3. Future of Heavy Flavor in STAR 

The Time of Flight system significantly enhances STAR's particle identification capa- 
bilities and with the large dataset recently taken we can improve all the measurements 
presented here. Further into the future, we plan to install a muon telescope system outside 
the magnet iron, supplementing the electron measurements with muons and resolving the 
T IS, 2S and 3S states. We are designing a heavy fiavor tracker consisting of a silicon 
pixel detector at 2.5 and 8 cm from the collision point, a silicon strip pad detector at 
14 cm and the silicon strip detector at 22 cm, expected to be completed for run 14. We 
will be able to reconstruct D'^ and mesons, measuring Raa- We will also measure 
Ac, checking if more charm quarks end up in baryons in A-A collisions. Identifying J/ip 
originating from a secondary vertex we can measure the feed-down from B mesons to 
3/ip- The direct measurement of D mesons will also allow us to measure V2 for charm. 
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2.4. V2 of Heavy Particles 



We have measured V2 for the hght nuclei d, t and "^He produced in Au-Au colhsions 
at y^snn = 200 GeV [9]. The V2 values are consistent with expectations from nucleon 
coalescence. For the multiply strange particles 0, S and 51 wc also observe a V2 consistent 
with coalescence, this time from quarks. By fitting a cocktail mixture to the tt+tt" 
invariant mass distribution we have measured the pt dependence of the meson V2 for 
the first time [10]. At high pt, above 2 GcV/c, the p V2 agrees with Kg V2. There may be 
a disagreement for low pt. This would be interesting because if the p is directly produced 
as a qq pair it should follow a coalescence for two quarks while if it is a recombination 
of two pions it should follow a coalescence for four quarks. 



Jets are produced by the fragmentation of hard-scattered partons. The scattering hap- 
pens early in the collision, usually creating a pair of partons back-to-back in </> but boosted 
along the beam direction, dependent on the a;s of the beam partons. One can get some 
information about jets by looking at high-pt particles, but because of fragmentation the 
leading particle is not always a good indication of the jet energy or direction. Better is to 
reconstruct jets explicitly using a jet-finding algorithm or implicitly through a technique 
such as two-particle correlations. Jet-finding algorithms group particles that came from 
fragmentation of a hard-scattered parton [11] [12] [13]. They have some assumptions built 
in though, such as the shape of the fragmentation pattern, and if jets are sufficiently mod- 
ified in a medium these reconstriiction algorithms could bias the measured jet properties. 
In A-A collisions there are typically many particles within the 'cone' radius that are un- 
related to the hard-scattered parton which must be subtracted. Two-particle correlations 
can include all pairs of fragments and the un-related background should be subtracted 
naturally. But there is some concern that other types of correlations can contribute. In 
the next sections we will describe observations of jets via two-particle autocorrelations 
followed by jets via explicit jet reconstruction algorithms. 

3.1. Two-Particle Correlations 

Two particle correlations measured by statistical measure 



are independent of 0i -|- (/>2 because of the 27r acceptance and nearly independent of 
rji + r]2 within the acceptance of the STAR TPC, so we can average over those variables 
without loss of information. The process of averaging over a stationary variable is called 
an autocorrelation. Averaging over 0i -|- (j)2 and 771 -|- r]2 results in a joint autocorrelation 
on {VAj^I^a), where rj^ = rji — 772 and (f)^ = 4>i ^ 4>2- To study angular correlations we 
can integrate over the entire observed pti and pt2 or a sub-range. Examples of joint 
autocorrelations for selected centralities from Au-Au collisions at ^snn = 200 GeV are 
shown in Fig. 2. We see evolution with centrality from a correlation structure consistent 
with p-p collisions for the most peripheral collisions to a structure highly elongated on 



3. Jets 



P(Pl,P2) - P(Pl)p(P2) 
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Fig. 2. Samples of {'ri/^,(f)^) joint autocorrelations for 200 GeV Au-Au collisions. The left-most panel is 
most peripheral (84-94%) while the right-most panel (5-9%) is nearly central. The middle two panels are 
on either side of the sharp transition, 55-64% and 46-55% centralities. 
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Fig. 3. We have subtracted the cos((/)a — i"): cos(2</>a) and 2D exponential components (determined by 
fitting) from the data leaving the same-side 2D Gaussian. The exponential affects only a few bins near 
(0,0). The away-side is flat, so we plot only the same-side, with 0^ and rj^ in a one-to-one aspect ratio. 
The circles are radius 0.7, indicating a typical jet-finding algorithm radius. 

TjA for central collisions. The peripheral correlation shape is well fit with a sharp 2D 
exponential (for HBT and e+e" afi'ecting only a few bins near (0,0)), a 2D Gaussian, an 
away side cos(<^a "TI") term, a ID Gaussian on rj^ (which is only important for peripheral 
centralities) and a constant offset. For mid centralities we also need a cos(2(/)a) term. 

We emphasize the centrality dependence of the 2D Gaussian peak. We fit each central- 
ity and subtract the 2D exponential, cos((/)a) and cos(20a) components from the data. 
The result is shown in Fig. 3. We see that the aspect ratio of the same-side peak changes 
quite dramatically, from a two-to-one ratio along the 0a direction in p-p collisions to a 
three-to-one ratio along the 77A direction in central Au-Au collisions. At the centrality 
where the aspect ratio changes from the (f)/^ direction to the 77A direction the 2D Gaussian 
amplitude increases dramatically as does the amplitude of the away-side cos((/)a — ti")- We 
interpret the same-side peak as the correlation of intra-jet fragments from hard-scattered 
partons and the away side cos(0a — t^) as the correlation of the inter-jet fragments. 

To support our interpretation of the same-side peak arising from hard scattering we 
examine its pt structure by selecting on pt of each pair included in the correlation. We 
find that for pairs with each particle around pt = 1 GeV/c the elongation along ryA 
happens at the same centrality and is just as strong as for all pairs. Having pairs of 
particles around 1 GeV/c from a thermal source should be very rare. 

Besides the study of jet-like structure we have also characterized the cos(2(/ia) compo- 
nent, {Ap[2]/ y/p^), that we extract from fits to 2D angular correlations. Previously we 
observed that the Ap[2\/ ^pj-d centrality dependence appears to be independent of beam 
energy and when plotted on relative impact parameter b/bo has a simple shape [14]. The 
beam energy and impact parameter dependence can be written as 

p{h)vl{2B}{h) = Ap[2]/^{ = 0.0045R(Vsr^)e2nbi„(b). 



p(b) is the particle angular density, nbin(b) is the number of N-N binary collisions (which 
depends on impact parameter) and we define u|{2D}(b). The important result here is 
that factor R{y/SNN) = k)g"200/i3 s?'^^'' depends only on beam energy and eopt(b) 
(the optical Glauber eccentricity) depends only on the initial geometry of the collision. 
We extend this to include the pt dependence of ?,'|{2D}(b) [15]. We measure two-particle 
correlations as a function of the pt of one of the particles. Extracting the cos(2(/)a) 
component we expect it should have a leading pt dependence which we remove. Defining 

Q(yt,b) = po(yt,b)^|{2D}(yt,b)/pt 
Po(b)t;|{2D}(b) (1/pt) 

so that Q(yt,b) is unit normalized we note that all centralities for Au-Au at ^snn = 
200 GeV fall on a single curve, defined as Qo(yt), which happens to be a Levy distribution. 
We can finally write the cos(20a) component in a form that is factorized on beam energy, 
centrality and pt- 

p(yt,b)?;|{2D}(yt,b) = {po(b)?;|{2D}(b)} {pt (1/pt) Qo(yt)} • 

3.2. Jets at Forward Angles 

In d-Au collisions at y^Snn = 200 GeV we can use a parton in the projectile deuteron to 
probe the gluon structure of the Au nucleus. Measuring at forward angles with 3 < < 4 
probes a range of x from 0.001 to 0.003. We measure the angular correlations of 7r° 
pairs. For peripheral d-Au we find prominent peaks on both the same-side and away- 
side, similar to what is observed in p-p collisions. For central d-Au collisions there is a 
strong suppression on the away side. The away-side peak behavior is consistent with the 
back-to-back jets in peripheral collisions but single jet production in central d-Au. The 
back-to-back jets are boosted to forward angles because of the parton x values. Single 
jet production is expected in color glass condensate theories [16], with a parton from the 
deuteron scattering off the low-x gluon field of the Au nucleus. The away-side suppression 
may also have explanations via other mechanisms [17]. 

3.3. Reconstructed Jets 

There is a variety of jet-finding algorithms which cluster particles within some 'cone' to 
measure properties of the parton whose fragmentation was observed. They don't use an 
actual cone, but they have a radius as a control parameter. We calibrate our jet finding 
using p-p collisions [18]. The cross section, shown in Fig. 4, is well described by pQCD all 
the way down to jet energies of 5 GeV. Jet fragmentation measurements are in reasonable 
agreement with Pythia plus Geant simulations. Using the away-side jet to measure kt in 
d-Au collisions we don't observe a modification of jet properties in cold nuclear matter. 

In central Au-Au collisions there will be many particles uncorrelated with the hard- 
scattered parton within any cone radius large enough to measure jet properties, and 
these particles need to be subtracted from the reconstructed jet. There are concerns that 
fluctuations in the uncorrelated particles can modify our reconstructed jet properties 
but we have studied this and believe we understand the uncertainties of the background 
subtraction. In Fig. 4 we show R;^^ and find that jets in Au-Au collisions are much less 
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Fig. 4. The left panel shows the jet cross section measured with STAR compared to a pQCD calculation 
for = 200 GeV p-p collisions. There is good agreement down to jet energies of 5 GeV. The right 
panel shows for jets reconstructed in central Au-Au collisions at -y/sNN = 200 GeV. Two different 

jet algorithms were used. The bands around the points represent the uncertainty due to background 
subtraction and the bar on the left is the uncertainty due to the energy scale of the reconstructed jet. 



suppressed than inclusive hadrons. Within the uncertainties there may be no suppres- 
sion. Jets are broader in central Au-Au than in p-p. We see this by running jet-finding 
algorithms with radius parameters of 0.2 and 0.4. In p-p collisions the larger radius jet 
contains about twice as much energy while in central Au-Au the larger radius jet contains 
four to five times as much energy. In two-particle correlations we saw that the same-side 
(intra-jct) component attributed to jets is significantly broadened along rj^ in central 
Au-Au collisions, in general agreement with this observation. 

We would like to use di-jets to study medium effects. There are possibly significant 
biases in our current understanding of di-jet reconstruction so for now we look at jet- 
hadron correlations. We determine the direction (and energy) of the reconstructed jet 
then look at 0jot — </>hadron for all hadrons within a selected pt range in the event. The 
signal to noise ratio is much larger for jet-hadron correlations than for combinatoric 
hadron-hadron correlations, so the details of the background subtraction are much less 
important. In Fig. 5 we show jet-hadron correlations for three associated hadron pt ranges 
and compare to p-p collisions. We see that the number of low-pt associated hadrons on 
the same side is larger for Au-Au than p-p while at high pt the numbers are equal. On 
the away side there are more associated hadrons at low pt in Au-Au while there are more 
associated hadrons at high pt in p-p. We quantify this in Fig. 6, showing that high-pt 
hadrons are suppressed in Au-Au relative to p-p while there is a significant enhancement 
of low-pt hadrons. This is similar to a two-component single-particle spectrum analysis 
taking the ratio of the hard components for Au-Au and p-p. The reason the inclusive 
Raa is small for low-pt hadrons is because Raa is dominated by soft particle production. 
Also shown in Fig. 6 is the energy transport. The away-side jet has less energy carried 
by high-pt particles, but most if not all of this energy loss is compensated by an increase 
of energy carried by lower-pt particles. We see that jets are substantially modified in the 
central Au-Au environment but not absorbed. 
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Fig. 5. Jet-hadron correlations for reconstructed jets with 20 < Pj^j.^;, < 40 GeV/c and three different 
associated hadron pt ranges. The solid symbols are for central Au-Au collisions at ^snn = 200 GeV 
while the open symbols for p-p collisions at = 200 GeV. This background subtraction assumed V2 = 0. 




Fig. 6. (Color online) The left panel is the away-side Iaa measured in central Au-Au collisions at 
Y^NN = 200 GeV. For all three reconstructed jet energies there is a suppression of associated-high-pt 
hadrons on the away-side (the highest pt jets having the least suppression for high pj^s>'°<:) and a sub- 
stantial enhancement of low-pt hadrons. The right panel, Daa = Pt{YAA ~ Ypp), shows that the loss 
of associated away-side energy in high-pt particles is largely compensated by an increase in energy in 
low-pt particles. 



4. Summary 

We presented measurements showing suppression of J/ijj (and maybe T) in Au-Au. 
Non-photonic electrons from D and B meson decays are also both suppressed. The in- 
stallation of the heavy flavor tracker will greatly improve our studies of quarkonia as 
well as allowing us to reconstruct D mesons to directly check for D suppression and flow. 
We will be able to reconstruct Ac and check if more charm goes into baryons in A-A 
collisions compared to p-p. 

We showed that at forward rapidities the away-side jet is suppressed in d-Au collisions. 
The 2D two-particle angular correlations show an evolution of the same-side jet peak 
from elongation in (f> for peripheral collisions to strong elongation in rj for central Au- 
Au collisions. Reconstructed jets show that in central Au-Au jets are not significantly 
suppressed but have fewer high-pt particles and many more low-pt particles. 
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